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ABSTRACT 


Deep water directional wave spectra, measured by a NDBC 
3-meter buoy off Monterey Bay, are transformed to shallow 
water using the linear refraction model by Dobson (1967). 
The transformed directional spectra are compared with 
measured spectra using pressure gauge arrays in shallow 
water at Marina and Santa Cruz. Refraction (Kr), shoaling 
(Ks) and Jacobian (J) transfer functions are computed. 

The modeling results and measured wave data at Marina 
are in good agreement for easterly waves at high frequen- 
cies. On the other hand, the linear refraction model gives 
highly variable values of wave solutions at Santa Cruz. It 
is believed that the unrealistic wave results at Santa Cruz 
are associated with the complex bathymetry of Monterey Bay. 
The present work examines the accuracy and limitations of 


linear refraction model by field observations. 
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ieee PRODUCT EON 


Wave climatology plays an important role in coastal 
dynamics as most energy is put into coastal region by wind- 
generated waves. Initially, ocean wave analysis was limited 
to a single wave train defined by a wave height and a period 
in a one-dimensional spectrum, but in the last couple of 
decades the emphasis has been placed on two-dimensional 
Spectra. 

The shoaling waves undergo refraction, reflection and 
diffraction mostly influenced by the wave approach angle and 
the bathymetry of the area. A spectral approach requires 
the superposition of a large number of progressive sinusoi- 
dal component waves. Therefore, waves of small amplitude 
(linear wave) are assumed. This assumption neglects 
reflection, diffraction and postulates long-crested waves. 

The area of Monterey Bay has been selected for this 
comparative study of ocean waves. This area presents an 
interesting problem for predicting directional wave spectra 
because of its variable shoreline orientation relative to 
prevalent eastward waves and the presence of the Monterey 
submarine canyon bisecting the bay. The Monterey submarine 
canyon is the largest in the western hemisphere and causes 


substantial perturbations on the nearshore wave field. 


Historically, Monterey Bay was one of the first, and 
presumably the most intensively studied, locations for wave 
rerract lone Wave refraction diagrams of Monterey Bay were 
first constructed in 1948 by Johnson, O'Brien and Isaac 
(Wiegel, 1964). The first intensive refraction studies were 
performed for amphibious landing exercises at Fort Ord 
(1954) followed by studies by the Corps of Engineers (1958) 
for small craft harbor at Santa Cruz. Wave refraction and 
littoral drift calculations have been performed by Dorman 
(1969) and Arnal et al. (1973). All these studies were 
performed using hand-drawn refraction diagrams. More 
recently, numerical refraction calculations have been made 
for a number of specific marine construction sites (e.g., 
Thornton (1978, 1980)) and erosion studies (Oradiwe, 1986; 
McGee, 1986). 

The objective of this study is to compare ocean wave 
measurements with the linear refraction model. The deep 
water two-dimensional spectrum is transformed to shallow 
water spectrum using transfer function computed from a 
linear wave refraction model by Dobson (1967). The deep 
water spectrum, So(f9,09) is determined from data acquired 
from the National Data Buoy Center (NDBC) using the Fourier 
analysis technique as proposed by Longuet-Higgins (1963) and 
Cartwright. ( Ieee The transformed spectrum S(f,%), is 


compared with measured shallow wave spectra for which the 


data had been acquired from pressure gauge arrays located at 
Marina and Santa Cruz. 

The background theory is described in Chapter II. The 
data acgquisition and equipment are mentioned in Chapter III. 


Chapters IV and V contain discussions and the conclusions. 


Il. THEORETICAL BACKGROUND 


A. WAVE REFRACTION 
1. General 

The refraction of ocean surface waves on the 
continental shelf is necessary for forecasting ocean waves 
in shallow water. Practical techniques for calculating 
refraction diagrams were developed by Johnson (1948) and 
Munk and Arthur (1952). The graphical technique developed 
by these investigators is called the Orthogonal Technique. 
With the advance of digital computers, the numerical methods 
were introduced by Griswold (1963), Wilson (1966) and Dobson 
(1967). The numerical calculations of wave orthogonals (or 
rays) can be quickly and accurately plotted for engineering 
purposes. The study of wave transformation that includes 
refraction over complex underwater topography and in large 
areaS 1S more commonly accomplished through numerical 
models. (Hayes, 1980) 

2. Refraction Theory 

In this study, linear wave theory is applied to 
shoaling and refracting waves, which neglects reflection and 
diffraction assuming energy propagates only in the direction 
in which the waves travel. It also neglects amplitude 
dispersion. The wave refraction is based on generalization 


of Snell's law for phase speed variation and conservation of 


energy. The phase speed of the wave depends upon wave 
period (or frequency) and the depth of water at particular 
locations. Wave energy is conserved between wave orthog- 
onals, where orthogonal is perpendicular to the wave crest. 
The orthogonal indicates the direction in which shoaling 
wave progresses. Based on linear wave theory, it is assumed 
that the wave period will remain constant but’ the 
wavelength, the velocity and the height, will vary as the 
wave shoals toward the beach. 

In general, the phase speed of a small-amplitude 


linear water wave is expressed as 


C2 = 2 tanh (kd) (2a) 
where g is acceleration of gravity, k the wave number and a 
is the depth of water at a particular location. Consequent- 
ly, in shallow water the phase speed of the wave varies with 
depth only and Snell's law becomes applicable in the case of 


straight and parallel contours as 


See. GC. (ZZ) 


where ois an angle between the wave crest and local depth 
contours and C is phase speed. The subscript zero refers to 


deep water conditions. 


The energy of a wave when it reaches the beach is a 
function of the wave height, and the wave height at the 
shore is related to the wave height in deep water. The 
variation of wave height within two wave rays is described 


(Mei, 1980) as 


Coie 2 20 2 (238 
pee Nite! bai! 
or 
H = Ho Kc Kr (2.4) 
where 
Cg 
Ko = lage = Shoaline Coetrreremes (2. 59 
b 
Kaa [es = "Refraction Cociricieme (2 7Gy 


where b is the distance between two rays and Cg iS group 


velocity and may also be represented in the form of 


2 kd 


SE Dal 


Cg = £ [1 + 


Using the relationship C = Cp tanh(kd) (Kinsman, 1965), 


Equation (2.5) may be wr@ttenmac 


k 1/2 2 kd ,-1/2 
x) «(+ Sinn 2kd! a! 


From Equation (2.2), the change in angle determines the 
change in wave crest or b value. Thus the value of Kp may 


be computed from 


a 
COS ap cos a4 
or 
Kr = 2041/2 = pa 9) 172 (2.6a) 
b cos «a 
where 
Sin a 
= Se) lic = ts 
0 


In any linear wave refraction model, shoaling and 
refraction coefficients are first computed and then the wave 
height is determined. 

3. Linear Wave Refraction Model 

Munk and Arthur (1952) described the theory of wave 
intensity along a refracted wave. Dobson (1967) followed 
the theory and developed a linear refraction model. The 
Dobson model is briefly described here. 

The definition sketch for the wave ray equation is 


Ewmowm In Figure 2.1 im that A,B, is tangential to circular 





Wave Rays 





Figure 2.1 Definition Sketch for Wave Ray 
Equations (Dobson, 1967) 


arcs of radius Cét and (c+é6C)ét and S and n are very small. 


Two basic characteristic equations for the propagation of a 


wave ray are the ray curvature equation 


and an equation for the ray separation factor 


1967) 














(2.8) 


8B (Dobson, 


me : Z 
ge + SB RS + ec HG = 0 C259") 
at Dn 
Dobson manipulated these equations using Gauss 
integral theory into a convenient form where Eqs. (2.8) and 
(2.9) can be easily applied in terms of cartesian 
coordinates and time variable as 
meee. % ace Be 
R = q@isin O ae Cos @ ay (2.10) 
a* de 
Bp a gg = 0 (25 t) 
at 
where 
O he 
Pee a2 Cos. o ae + sina =~] (222 } 
OX oy 
2 D Z 
’ 0 ; fe cee: 
oon) = C[sin"s —§ - 2 sin « cos a = + cosa 5) (22135) 
3x i oY 


The phase speed is a function of the coordinates so 
thats C [= “Gee. )] The wave ray can be constructed with 
Equation (2.10) within a coordinate system. The ray curva- 
ture R is computed using finite difference methods between 


two consecutive points along a wave ray by iterating 


An = RAS (2.8 


at each step. So the new coordinates can be computed as 


OW oe OS es a (2. ee 
Yours = Yo + AS cos @ (2 ey 
Whee: 
= seers 
or 2 
—- _ % * “g+ag 
—— 2 


The ray path is constructed grid to grid over the sea bottom 
until it reaches the required location. 

The ray separation factor 8 in Equation (2.171). 
computed applying Fox's formula for a second order linear 
equation. This is a finite difference method and requires 


that the separation factor be Known at the starting and 


10 


preceding points. This condition is conveniently satisfied 
by starting the ray in the deep water where ¢ is equal to 1. 

The characteristics of the wave at any point in the 
coordinate system may be determined from the solution of 
Equations (2.10) and (2.11) provided the speed, C(x,y) is 
known. However it is common to consider waves of different 
deep water characteristics travelling towards the coastline 
over a known bathymetry, so that it is more convenient to 
consider the problem in terms of water depth, d(x,y) rather 
than wave speed. Thus these equations are expressed as 
functions of the mean water depth. 

Dobson used a quadratic least square surface to fit 


depth locally at intermediate points using 


Gd = 04 + Pox + pay + Pyx? + Poxy + Sey? (27) 


where the partial derivative of depth, da, with respect to x 
and y is expressed as functions of 6. Therefore the local 
water depth can be obtained by interpolation method. 

The bathymetry used in this refraction program was 
obtained from the original NOAA data and projected onto a 
six second (latitude and longitude grid) modified Universal 


Transverse Mercator (UTM) grid after screening the data for 


errors and bad points. Insufficient data were available for 
off-shore bathymetry. Thus it required interpolation to 
fill the empty grid cells. Piece-wise linear triangulation 


all 


was used to interpolate to a 200 meter rectangular grid. A 
nine-point weighted linear average was applied to smooth 


bathymetry (Oradiwe, 1986; McGee, 1986). 


E.. DIRECTIONAL WAVE SPE@GReur 
1. General 
Directional characteristics of ocean waves are 
important, not only for understanding wave generation 
mechanisms but also for practical problems such as wave 


forecasts in the design of coast and harbor protection 


works. Directionality 1s equally important in the 
refraction of ocean waves. Various techniques incorporating 
floating buoys, pressure gauge arrays, electromagnetic 


current meters and SAR have been utilized to determine 
directional properties of the ocean waves. In this stud 
data have been acguired from the NDBC heave/pitch/roll buoy 
for deep water spectra and pressure gauge arrays for shallow 
water. The theoretical background of these measurements is 
presented in the following sections. 
2. Deep Water Wave Spectrum 

Longuet-Higgins (1963) suggested a method to compute 
directional wave spectrum from water surface elevation and 
two orthogonal components of wave surface slope. The 
surface elevation is supposed to be the result of the 
Superposition of a large number of long-crested progressive 
sinusoidal component waves in random phase (Battjes, 1972). 


Several other investigators (e.g., Cartwright, 1963; Long, 


eZ 


1980; Mitsuyasu, 1975; and Long and Hasselmann, 1979) have 
suggested and refined the technique. 
The directional spectrum S(f,®) as given by Longuet- 


Higgins (1963) is expressed in the form 
S(f,0%) = S(f) ° D(f,a) (28) 


where S(f) is the one-dimensional frequency spectrum and 
mer, o) 1s angular distribution function with the condition 
27 


f D(f,a)do = 1 (2.19) 
0 


Equation (2.19) can be expanded by a Fourier series 


D(£,2) =2 


cae ) {a cos(no) + b sin(nx) }] (252 00) 

2 n n 
n=l 

The coefficients of the first two harmonics can be estimated 

from the co- and quad-spectra from heave (index 1), and 

slope signals (indices (1,2), correspond to x,y slopes) 


(Longuet-Higgins, 1963; Long, 1980) as 


27 Q19 
a = (Dac) cos o dou= x (ew) 
0 
Pee D(t,e)esin« daw = 722) 
0 


i Ie 


27 Cc - C 





ll = | D.Gh en) COSN 2 Gai oa a a (2.228 
2 0 K*C 
cee 
bo DUE, 0) samecmeclamee = (2.225 
0 K Cia 


where the Constante ss Genteea as 


SO)" Obs 


alg 


41/2 


A heave/pitch/roll buoy and pressure sensor array 
provide a measure © 7 a0. aude Ny at fixed points. The co- 
and quad-spectra of pairs between these three quantities are 
computed using the Fast Fourier Transform. 

The wave directionalities are determined by the 
directional distribution function. Equation (2.20) may be 


rewritten as 


[r,Cos(a-a,) + r,Cos 2(a-a5) J (2.259 


where 


In the case of unimodal spectrum, Longuet-Higgins 


(1963) proposed a cosine powered distribution function 
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icin) mG sme COS = ) (2526) 
where 
© 426-1 S41) _ 425-1 (8!)" 
S(2)) - Tl 2S (2S) ! 
is a normalizing function to satisfy Equation (2.19). The 


parameter S is a function of wave frequency and « is the 


mean direction computed from 


Cartwright (1963) showed that the value of S is determined 





by 
qi 
il 
S — 
= (1 - r,) 
L cayas aes (x4 + l4r, + i aac 
S32 = = 
Zl r.) 
The cosine distribution model is used most in uni- 
modal wave analysis. The same expression is used in this 
study. 


3. Predicted Shallow Water Wave Spectrum 


The shallow water wave spectrum can be predicted 


from the deep water spectrum by multiplying the squares of 


i 


both refraction (Kp) and shoaling (Kg) coefficients and the 
Jacobian of the direction function by each frequency 
component. Hence the deep water spectrum is transformed 


into shallow water spectrum by 
D2 
2 (roe) = So (fo,%o9)Kp Kg J (2.299) 


where the Jacobian 


£(£ 51%) 
J a ACE os 


Since the wave period remains invariant, J only depends upon 
the changes of local wave direction. 
4. Measured Shallow Water Wave Spectrum 

The two-dimensional shallow water wave spectrum is 
determined from the data acquired from shallow water 
sensors. The same technique as described previously for 
deep water spectrum 1s used. In this study, the pressure 
gauge arrays have been utilized in which subsurface pressure 
fluctuations are measured to infer the surface elevation 
(Kim, 1974; Pawka et al., 1976). 

From linear wave theory, the relation between wave 


elevation (") and subsurface pressure (P) 1s described by 


cosh kd 


n(t) = Laan oe Pie) (2.268 
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where 


ad = water depth in meters, 
Z = vertical distance from mean sea level to 
pressure gauge (Z is positive upward from 
sea level), 
P = pressure measured in meters of water head 
about a zero-mean pressure. 
The term in the brackets in Equation (2.28) is the spectral 
transfer function relating pressure to surface elevation as 
a function of frequency and depth. Equation (2.28) can be 
modified to include differential distance and elevation 


between two sensors of an array to compute co- and quad- 


spectra of surface slope as 


cosh kd 2 


.~ 1 ,_cosh kd 2429 
Za, = > [eee emp one! ( 
oe AX BA 
ax 

Jk cosh kd ‘ (ood) 
fe eel eree(aan)) op, Ape 
dx Ti 


Bere OPS is subsurface pressure adifference between two 
sensors on the x-axis and Ax is differential distance along 
the same axis. Similarly, pressure power spectra and co and 
quad-spectra can be computed along y-axis. Once the power 
spectra, co, and quad-spectra are obtained, the directional 
spectral density may be determined using the same technique 


outlined earlier for deep water directional spectrum. 


iy 


Iitl. DATA ACQUISITION AND SENSORS 


Age LOCATION AND TIME 

The area of Monterey Bay was selected as the study area 
which poses complex problems due to the highly irregular 
bathymetry of the Monterey submarine canyon. Three sensors 
were utilized for acquisition of wave data. The locations 
of the sensors are as follows (see Figure 3.1 for bathymetry 


and Location 


Noe Type of Sensor Location Depth (M) 

(ey NDBC Fol l/pire 36°48'00''N 2003 
wave buoy 1229°23'59' 'W 

Dis Array of pressure 36P°42.0'N 15.0 
gauges (4) 
Marina 121°48.9'W 

oe Array of pressure 36°57.0'N T Seale 
gauges (4) 
Santa Cruz 122°00.2'W 


Selected wave cases from January 1988 are used in the 
analysis. During the winter season, weather conditions are 
generally favorable for creating higher energy waves. In 
addition, an abnormal storm was encountered in the region. 
To select a particular date and time, the conditions for 
narrow band wave spectrum were emphasized to simplify the 


interpretations. 
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In comparing coincident data sets, on the average a time 
lag of one hour was assumed for the waves to travel from the 


deep water buoy to the shallow water array sites. 


ion DEEP WATER WAVE DATA 

Deep water wave data for the pitch/roll buoy were 
obtained from the NDBC. The NDBC 3-meter buoy measures 
heave, pitch and roll movements with respect to 3-axes 
magnetometer. The magnetometer axes are aligned with fore 
and after (bow positive), athwartship (starboard positive) 
and vertical (positive upward) axes of a buoy. The analog 
sensor measures vertical acceleration and displacement in 
addition to hull pitch and roll. The magnetometer bow and 
starboard output, together with vertical displacement, pitch 
and roll output from analog sensors are sent to the 
Directional Wave Data Analyzer (DWDA) after voltage 
amplification and filtering. 

The hull azimuth angle is calculated from bow and star- 
board magnetic components with necessary corrections made 
for pitch, roll and hull magnetic fields. The magnetic azi- 
muth is converted into true azimuth using the magnetic vari- 
ation for the site. The true azimuth is then used to 
calculate East-West and North-South components of hull 
slope. One hertz time record of hull displacement and 
Slopes are stored in DWDA for spectral analysis and trans- 


mission to shore station through GOES (Steele et al., 1985). 


20 


Data provided on magnetic tape included environmental 
parameters, wave spectra, directional wave parameters and co 
and quad spectra. Directional wave parameter data were 
utilized to determine directional wave spectrum. The 
unimodality of the directional waves was confirmed by 
studying the regional weather maps to conveniently apply 
cosine angle distribution function. Two-dimensional spectra 


are shown in Figure 3.2. 


C. SHALLOW WATER WAVE DATA 

The calibrated data for arrays of pressure gauges were 
obtained from the Scripps Institute of Oceanography. The 
array consists of four pressure gauges located in a six- 
meter square area which measures pressure in centimeters of 
water. Each sensor is located about 50 to 100 cm above the 
bottom. Pressure power spectrum is calculated which in turn 
1s transformed into surface elevation and surface slopes 
ny and ny, spectra by applying linear wave theory transfer 
function as described in Chapter II. 

For the computation of the cross-spectra, the following 


mamues were utilized: 


Ct =S2Sec, 


sampling interval 


total data points = 2048, 


- data points for one sample = 128, 
Beecocal record length = 32 min, 


- Nyquist frequency = 0.25 Hz 
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Part of Two-dimensional Deep Water Wave 
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Ze 


- frequency bandwidth = Af = .00391 Hz, 
- degrees of freedom = v = 32. 

The spectra were cut at frequency of 0.2 Hz in the 
higher frequency range due to spurious overestimation of the 
mranstere function. The wave period for high quality 
estimate is approximately 5-20 seconds. 

To ensure accuracy and to check computations, the 
results were compared and verified with the published data 
in the Coastal Data Information Program of January, 1988. 
Two-dimensional spectrum at the peak energy frequency is 
depicted in Figure 3.3 for the Marina array on 29 January 


Boia : 


D. OCEAN WAVE TRANSFORMATION PARAMETERS 

The Dobson linear refraction model has been used to 
estimate the shoaling transformation of the wave direction- 
mieespectrum. The reverse projection method is employed. 
That is, the refraction program was run backward from 
shallow to deep water, commencing at a shallow water loca- 
tion in required depth, for Marina 15 meters and for Santa 
Cruz 13.1 meters. Rays of a particular frequency were 
propagated offshore in increments of .2 degrees (0.2°) over 
the range of all possible incoming wave angles. Once the 
depth is deep, the rays are stopped and deep water wave 
directions are recorded. The rays were then returned to 
retrace the same path to verify the backward refraction and 


memealculate transfer coefficients Kp, Kg and J. 
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Figure 3.3 Two-dimensional Shallow Water Wave 
Spectrum Measured at Marina on 29 
January 1988 
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Lome nelude = rehraemrenm, Over the shelf to the north of 
Santa Cruz, supplemental bathymetric data were meshed into 
the model. The shelf was more coarsely gridded to a 1 Km 
grid because of the large area involved and due to the fact 
that the waves are primarily in deeper water while 
traversing the shelf. The transfer parameters calculated 
from the two areas were multiplied together to determine a 
close relationship for use in spectral transformation. Some 
example ray traces are shown in Figure 3.4. 

1. Shoaling Coefficient 


Shoaling coefficients, Ke computed from the model 


for Marina and Santa Cruz sites are depicted in Table 3.1. 


TABLE 3.1 
SHOALING COEFFICIENTS, Ke 
Freq. Shoaling Coefficients 
(Hz) Marina Santa Cruz 
(15 m) (1 om) 
0.030 Zee Oro ok 
0.040 1.689 1.796 
0.050 Nip 63 )8, ee 7 2 
0.060 eo] i263 
O20 70 1.068 fl See 
0.080 0.979 i022 
0.090 0.918 0.951 
0.100 O37 7 0.902 
Os 10 Ores a4. 0.868 
0.120 O7os 7 O23847 
O23 0 O-se4 Oma36 
0.140 0.839 0.834 
0.150 0.850 0.839 
0.160 0.867 Os 50 
O2 170 Ooo 7 O58 6 5 
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Figure 3.4 Ray Traces from Marina to Deep Water 
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2. Refraction Cociiterene 

The model computes refraction coefficient, Kp, seem 
every 0.2° of nearshore angle. Kp values were averaged jiam 
every degree of deep water direction. In the linear refrac- 
tion model, an assumption is made that the underwater topog- 
raphy changes are small and gradual but, occasionally, this 
assumption does not hold in Monterey Bay because of the 
canyon. As a result, the refraction model sometimes gives 
erratic estimates of “the @kp values. Therefore, all i@ 
values above three were truncated to the value of three. [It 
is also observed that refraction estimates are more stable 
in higher frequencies (short period waves) than lower 
frequencies (long period waves). See Figure 3.5. 

36 .Jacobuan lita) ; 

The Jacobian (a is calculated from the change in 
the deep water direction as a function of near shore angles 
for each wave frequency for all possible wave ray direc- 
tions. A plot of og versus % (Figure 3.6) shows a slope of 


nearly 45° over most of the angles close to the beach 


Normal. The slope tends to increase as the wave ray 
undergoes greater refraction. The Jacobian calculated from 
these values (Figure 3.6) gives Similar results, 1.e., the 


Jacobian estimates are constant and closer to one when the 
approach angle is within 25°-309° of the beach normal and 
increasingly unstable/erratic at extreme limits. Values at 


the limits may not be useful for computation of directional 
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spectra particularly in the lower frequency range of .07 Hz 
and less. 

The Jacobian appears to behave similar to the 
refraction coefficients. The Jacobian estimates at higher 
frequencies are more stable than at lower frequencies. 
Figure 3.7 shows the Jacobian estimates at different fre- 
quencies for the Marina array location. 

4. Transformation of Deep Wave Spectrum 

The deep wave spectrum is transformed to shallow 
water wave spectrum using Equation (2.27). Transformation 
parameters calculated from the refraction model are 
multiplied with deep water directional spectrum to obtain 
shallow wave directional spectrum. Figure 3.8 depicts the 
transformation of 0.09 Hz wave for different wave 
directions. The graph of Kcg is not shown since at a 
particular frequency it has a constant value for all 


aeerections. 
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Figure 3.8 Transformation of Spectrum from Deep Water 
to Shallow Water at Marina for Peak 
Frequency 0.09 Hz 
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IV. RESULTS AND DISCUSSION 


A. REFRACTION MODEL PREDICTIONS 

Transformation parameters, Kp and Jacobian, determined 
by the linear wave refraction model were variable, with a 
tendency to increased variability with increased refraction. 


The lower frequency waves are more refracted than high 


frequency waves. Therefore, these coefficients tend to vary 
smoothly in higher frequency (short wave period) and 
increased in variation with a decrease in frequency. The 


variability of Kp and J is greater for the Santa Cruz site 
for which incident waves must undergo significant refrac- 
tion (see Figures 4.1 and 4.2). The characteristics of each 
location and refraction model predictions are described 
here. 

The Marina location is shadowed by Point Pino in the 
Southwest and Point Santa Cruz in the Northwest. High fre- 
quency waves, 0.11 Hz and above, are limited to directions 
from 2359-3109 at the array site, whereas low frequenes 
waveS May arrive from even southerly directions after 
considerable refraction (refer to Figure 3.4). Thus, the 
overall cut-off limits are 235° in the South and 310° ites 
North. The wave rays, at low frequency, are sparse and tend 
to form caustics which result in higher variability for 


estimates of coefficients. The variation of J also depends 
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at Marina and Santa Cruz 
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upon the direction of approaching waves in addition to the 
frequency. If the wave ray coincides with, or is closer to 
the beach normal, the J values smoothly vary over the 
directions. As the angle of approach increases, J tends to 
vary with increasingly high estimates. 

The Santa Cruz site poses more complex problems than 
Marina because of its location. The waves approaching from 
the westerly direction makes almost a 90° turn to arrive at 
the shoreline. The low frequency waves travelling north and 
north-eastward may arrive at this location. The ray pattern 
shows (Figure 4.3 a,b) scattered and crossing rays (caus- 
tics) with occasional or a few wave rays from the west or 
higher approach angle at low frequency. This gives very 
large and highly variable estimates of J. The ray pattern 
becomes more regular with increase in frequency. It is also 
noticeable from the ray patterns that the J estimates are 
generally large and more variable compared with the Marina 
values (Figure 4.2). The results of Kp estimates behave 


Samilarly. 


Ee WAVE DATA SELECTION 

Data for wave measurements for the month of January 1988 
were examined. A number of cases were analyzed out of which 
three are described here. The selected days are 4 January 
Peewee rPstT, 16 January at 0200 PST and 29 January at 2000 PST. 
These periods provide narrow band energy density at differ- 


ent frequencies and direction. Figures 4.4, 4.5, 4.6 and 


3, 
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(a) Low Frequency 


Figure 4.3 Wave Ray Traces at Santa Cruz 
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Tables 4.1 and 4.2 


provide a 


comparison 


of measured 


directional spectra at the three locations for each day. 


TABLE 4.1 


DEEP WATER ENERGY DENSITY WITH MEAN DIRECTIONS 


FREQ 3 JAN O800 PST 18 JAN 0200 PST 29 JAN 2000 PST 

(HZ) DIR ENERGY DIR ENERGY Ne) Jonge ENERGY 
(1/sec) (Deg) (m2-sec) (Deg) (m2-sec) (Deg) (m2-sec) 
0.0300 17324 0.000 283.4 0.000 S373 0.000 
0.0400 Goes 0.000 297.3 Oo. 000 aoa. 3 0.000 
0.0500 240.5 0.000 265). 0 0.000 PAPE 0.000 
0.0600 e214. 2 0.000 ewe. 5 Pomoc 224.1 0.804 
wee? OO 2G. 6 24) Zo. 6 V2. 7 OO POG <5 PAB alliPabe. 
0.0800 270e0 ie a69 292 80 pO. 196 Pe See 6.346 
0.0900 ESS lee eke Ole «2 Cee onen | 2 die Se ore 
0.1000 242.8 io 10 S lieve Ps oI ii 20S 44256 
o- 1100 23O7 2.398 303.7 iyo 36 27.656 3.007 
oO. 1200 223.3 2755 i eo dia Wes) 2, 2 2enk 2,329 
3 OO eat oad, 2vo5l Bees 10.475 POS mas de 7/8, 
0.1400 a2 3.9 2.190 3136 8.089 269.5 Ji zeleieS 
0.1500 2 FG 1.699 6) ee | eres 2s 265.7 Ors6.5 
0.1600 2 iOe 1 1.418 Bee. oS PHA S28 ZOO MO 0.609 
0.1700 row. 3 1.347 Sree 2315 259.3 0.350 
O18 00 ef Oi seZ 0.889 309.6 1.940 25976 Once 
0.1900 i? D5.0 0.899 3:7 ee if 5.5.0 26 faae 0.184 
0.2000 ler eS mms Oo 34 oie ae Ia, Bao se 0.146 
te 100 Gee 0.670 Sale 8 0707 2 355.6 0.166 

TABLE 4.2 
WAVE CHARACTERISTICS OF SELECTED CASES 

SacACTERISTICS 3 JAN 88 18 JAN 88 29 JAN 88 
Hsign. epee. Ii Jc) ju 2 
FreQneak O77 1er iz 0.07 Hz UaO2 HZ 
FreQmean Ons Hz 0.10 Hz Gei2 H2 
Pigeectionnzean 220- 5010" DUS es 
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The California coast was hit by an unusually severe 


storm during 16-18 January 1988 which resulted in high 


waves, strong winds and sea level well in excess of 
predicted values (Cayan, et al., 1988). Widespread damage 
occurred along the Southern California coast. This storm 


developed about 300 NM west of San Francisco on 16 January 
1988 and traveled south-east at abou 33 5hnocc. The storm 
center passed off Monterey Bay on the morning of 17 January 
1988 continuing on its SE track and making landfall at Avila 
Beach in Central California. The bay experienced record 
high waves 14 hours after the storm passed from the area. A 
Significant wave height of 7.9 meters was recorded on the 
morning of 18 January 1988. The dominant wave period was 
about 14 sec with a direction of 284° with narrow =a 
energy measured at deep water and shallow water locations 
(Fagqure 4. 4)¢ 

Deep water energy density with mean direction for 
selected cases are given in Tables 4.1 and 4.2. On 29 
January 1988 measured spectral density of deep water waves 
and shallow coastal site are shown in Figure 4.5. 

Reviewing Tables 4.1 and 4.2 and other published data, 
deep water waves generally travel eastward CGuring the month 
of January with slight variability of about 20° on either 
side of true East. This direction is dominant in the 
energetic middle frequency range of 0.08 to 0.13 Hz. At 


other frequencies, the direction tends to be variable. 


Gx 


The third case of 3 January 1988 provides a different 
view from the normal as it contains sea wave energy at 0.13 
Hz frequency with a mean southerly direction of 228° (see 
mereure 4.6). In this case, deep water waves undergo less 


refraction to arrive at the Santa Cruz array. 


C. COMPARISON OF MEASURED AND CALCULATED SPECTRA 

The directional spectra from wave measurement data were 
determined using a cosine-powered distribution function 
feguation (2.26)). The value of the S parameter, in most 
cases, was less than 20 and relative directional energy was 
in a broad band of 50° to 60° spread. 

Predicted and measured shallow wave directional spectra 
for Marina are compared in Figures 4.7-4.9. The deep water 
wave approaching directions are limited between 235° to 
PaO? . These limits change slightly with decrease in wave 
frequency. In the low frequency range, the predicted 
spectral energy 1S overestimated (see Figures 4.7 and 4.8). 
The deep water spectra of high frequency waves from the West 
can be transformed into shallow water with reasonably good 
estimation comparable with actual wave measurements. The 
predicted spectrum from westerly waves is shown in Figure 
4.9. The predicted mean shallow water direction tended to 
be more aligned with the deep water spectrum than the 
measured shallow water spectrum. 

The transformation does not appear to appropriately 


shift the mean direction of spectral peak energy of deep 
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water wave spectra. Two additional cases of 5 January 1988 
and 30 January 1988 are shown in Figures 4.10 and 4.11 for 
high frequency spectra. The mean approach direction for 


these frequencies in deep water and shallow water are: 


Deep Shallow 
5 January 1988 0.12 Hz 255° 272 
30 January 1988 0.13 Hz 292° 281° 


The predicted spectrum in shallow water also indicates high 
Spectral density in the mean direction. 

At higher frequency sea waves, the refraction model 
overpredicts spectral density by less than 18% when summed 
over all angles at a particular frequency. For swell waves, 
the results are largely exaggerated. 


Results of the Santa Cruz location are depicted in 


Figures 4.12 and 4.13. The results are quite different from 
Marina. The predicted spectra are highly overestimated at 
lower as well as higher frequencies. The relatively simple 


case of 3 January 1988 was particularly selected for Santa 
Cruz to consider ocean waves from the southwest direction. 
Even this simple case does not provide satisfactory results. 

These differences are broadly associated with bathy- 
metry. The non-linear effects due to diffraction and 
currents were neglected which may introduce errors. It is 
surmised that the complexity of bathymetry in Monterey Bay 
may contribute significantly to overestimation and variation 


of transformation parameters. 
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Deep water directional wave spectrum was determined from 
data acquired from the NDBC buoy and transformed to shallow 
water using the linear wave refraction model. The 
transformed spectrum was compared with measured shallow wave 
spectra at Marina and Santa Cruz. 

The transfer parameters of refraction coefficients and 
Jacobian predicted by the linear refraction model are highly 
dependent upon wave period and direction of approaching 
waves. The Jacobian appears to be more sensitive to wave 
period and direction than the refraction coefficient. These 
coefficients are generally stable at higher frequencies 
and/or when the wave ray 1S normal or nearly normal to the 
coastline. 

The model tends to overestimate the transformed 
spectrum. The estimates are less than 18% in excess of the 
measured spectral density for easterly waves at frequencies 
of 0.09 Hz and higher for the Marina site. At lower 
frequencies the results are largely exaggerated. 

The model provides unrealistic large variations in wave 
height calculations for waves at Santa Cruz. This Ge 
presumably due to its location, orientation of shoreline 
with prevailing north-easterly wave direction and the 


presence of Monterey submarine canyon. 
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It is clear that the strong irregularities of bathymetry 
in Monterey Bay causes non-linear effects. The linear 
refraction model is not suitable for handling such problems. 


Further modeling efforts including non-linear effects are 


needed for Monterey Bay wave study. 
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